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Abstract 

We propose a test to compare the speed of muons and muon neutrinos produced 
in the same bundles in the atmosphere. By comparing the arrival times of primary 
muons and muons induced by neutrinos, experiments such as Super-Kamiokande, 
IceCube and ANTARES could verify recent hints on the possibility that neutrinos 
have a velocity larger than the speed of light. For smaller detectors such as MINOS 
and OPERA, where it may be challenging to recognize downgoing muons induced 
by neutrinos inside the active volume, another measurement may be possible thanks 
to the ns-time resolution. The arrival times of muons in a bundle produced by one 
cosmic ray interaction in the atmosphere should return a value compatible with 
the time resolution, while if neutrinos have larger velocity than muons some of the 
muons may arrive slightly earlier in a bundle of particles than others. 

Key words: 

PACS: Neutrinos, flux, kaons, pions 



1 Introduction 

With long baseline neutrino experiments it is possible to measure the neutrino 
velocity by the time difference between neutrino detection and production. 
Recently the OPERA neutrino experiment on the 730 km CERN-Gran Sasso 
neutrino beam has measured a difference of 60.7±6.9(stat) ±7A(sys) ns with 
respect to the time computed assuming that muon neutrinos are traveling 
with the speed of light observed by OPERA. There is no apparent energy 
dependence of this time difference between ~ 10 — 40 GeV, the energy spread 
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of the CNGS beam. This result, if confirmed, would imply that neutrino are 
tachyons with speed above v the speed of light by ' !t=£ = 2.48 ± 0. 28 (stat) ± 
0.30(sys) x 1(T 5 pQ. 

This result is in agreement with the result of the MINOS experiment that 
applied a similar technique but with a much larger error [2] . The OPERA result 
is however in disagreement with the observation of anti-electron neutrinos of 
energies of about 10 MeV from the famous supernova SN1987A providing 
an upper limit of of \v — c\/c < 2 x 1CT 9 |3j. If the effect is real, then this 
could hint to a possible energy dependent effect or other yet unknown flavor 
dependent effect. 

The interaction of a primary cosmic ray with the atmosphere produces a cas- 
cade with many kind of particles, and in particular neutrinos and muons. 
Muon neutrinos and muons are produced mainly via the decay of charged 
pions and kaons produced in the primary cosmic ray interactions and above 
about 10 TeV they can come also from prompt decays of charmed hadrons. 
This last component has not yet been observed. The neutrinos in this cas- 
cade, if the OPERA result is correct, should arrive before the muons from 
the same parent decay with a time delay that should change according to the 
neutrino path length that depends on its zenith angle 9. In underground de- 
tectors muon neutrinos are detected looking for induced muons produced by 
neutrino charged current interactions in the rock, and the ice around or inside 
the instrumented region. Hence, a time spread should be observed between the 
muons produced by the pion or kaon decay directly and the muon produced 
by neutrino interactions. 

The path length from the meson decay point is a few tens of kilometers for 
vertical neutrinos and up to ~ 300 km for near horizontal neutrinos. Assuming 
that the time difference observed in OPERA is not energy dependent and it is 
60 ns, nearly horizontal neutrinos should arrive up to 28 nsec before the other 
secondaries. In Ref. [3j a table of average production heights neutrinos in the 
atmosphere has been reported. The typical production height for neutrinos 
of energy above 20 GeV can be 17.6 km at the vertical, 94.9 km at cos# = 
0.25 and 335.7 km at cos# = 0.05, which would correspond to 1.4, 78 and 
27.6 ns. The typical time resolutions of photomultipliers (PMTs), including 
the associated electronics, is between At ~ 1 — 3 ns [5f6"lP7) . Typical tracks are 
reconstructed using a number of PMT time hits TV and the time uncertainty 
associated to a muon track depends such number as At/yN. 

There are already limits to tachions in cosmic rays. The limits are obtained 
looking for example to signals before the main front of the electromagetic 
shower. But this kind of searches stopped some time ago and the last particle 
data book review related to this tachion searches is the one of 1994 [5]. The 
limits obtained are of small interest in the framework of the OPERA result. 
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However, if neutrinos are tachions, it is likely that other kind of tachions could 
exist and this search in very high energy cosmic rays could have a new interest. 



The MACRO experiment has done several searches for possible anomalies 
of the time differences between muons PITT] . The searches were dedicated 
mainly to study time differences of the order of a few msec or more, but 
Ref. [9] contains also the study of time differences at the ns-level. The results 
of this study was that the largest time differences between two muons was 
18 ns and that the time distribution was in agreement with the predictions. 
The statistics was limited to 35832 tracks in events with two or more tracks. 
In 1992 none was thinking to tachionic neutrinos and therefore there was no 
estimate of the number of tracks due to down-going neutrino together with a 
primary muons. In Ref [TO] this study was extended to about 140000 tracks of 
multi muon events, corresponding to about 4% of the total MACRO statistics. 
The time distribution was in agreement with the predictions. 



Currently the largest neutrino detector in operation is the IceCube detector 
at the South Pole. We will discuss in the following paragraphs the possibility 
to measures the time delay between the muon produced by the neutrino in- 
teraction in the ice and the muon produced directly in the atmosphere using 
IceCube and the compact detector inside IceCube called DeepCore. DeepCore, 
vetoed by the surrounding IceCube strings, has the possibility to separate 
events produced internally by neutrinos and the muons directly produced in 
the decays in the atmosphere. The real analysis is made difficult by the scarce 
density of photodetectors that prevents a reliable tracking of multiple muons in 
one event but the veto possibility offered by the PMTs surrounding DeepCore 
and the large dimensions of IceCube may make this measurement interest- 
ing. In fact the tracking can be applied to clusters of hits recognized by 'hit 
cleaning' algorithms that can be applied to recognize different topologies in the 
same time frame of the data acquisition. ANTARES has better time resolution 
for muon tracks than IceCube due to the smaller amount of light scattering. 
A similar analysis may be attempted there though the identification of the 
neutrino may be more tricky in the absence of a veto. Super-Kamiokande can 
measure internal neutrino events and has good pattern recognition capabilities 
as well as OPERA and MINOS that have also excellent time of flight mea- 
surement capabilities. Even if these detectors have limited dimension and may 
not be able to identify internal downgoing neutrinos, they have good track- 
ing capabilities and may perform a measurement of muon bundle timing. All 
muons should arrive at the same time compatibly with the time resolution, 
but muons induced by neutrinos could arrive earlier if they are tachions. 
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Table 1 

Parameters for the 5 degree polynomial in x = cos 6 and y = logioE as described 



by eq. [T] for the calculation in Ref. |13|12j for muon neutrinos and anti-neutrinos. 
2 The semi-analytical atmospheric neutrino flux calculation 



We have calculated the atmospheric muon neutrino flux using the tables pro- 
vided by the Japanese group in [12J. Their calculation is described in [13J. We 
have fit these tables for muon neutrinos and ant i- neutrinos with a 5 degree- 
polynomial function of (x,y) = (cos 9, log^E), where E and 9 are the energy 
and nadir angle of the neutrino respectively: 

(dN v \ 5 5 4 

l ° 9w [ ~dW I = ? Vix%l + Po + ? ViyV% + ? Pix v x%y5 1 W 



The parameters of the fit are given in Tab. [TJ 



Between energies of ~ 10 — 100 GeV the flux of atmospheric neutrinos is not 
a single power law but above about 500 GeV the flux can be described by the 
analytical function that we reproduce here from Ref. [T] 



dN Vll 
dlnE 



A,. 



+ 



1 + B P E cos 9* /e w 1 + B k E cos 9* /e k 



(2) 



where cos#* is equal to cos# for cos 9 < 0.3 and for horizontal directions 
(cos# > 0.3) has a larger value that accounts for the curvature of the Earth. 
Notice that we consider the neutrino fluxes at these energies up-down sym- 
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Table 2 

Parameters for the cos 0* function described in eq. [3j 
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Table 3 



Parameters for the high energy functions in eq. [2] for and fluxes EdN/dE in 
units cm -2 s _1 sr^ 1 , with the energy in TeV for HKKM2006 [IB] as derived fitting 
the tables in the energy region 500 GeV-100 TeV. 

metric. Above 500 GeV, muon decay does not contribute to the neutrino flux 
at these energies and scaling of hadron-nucleus interactions starts to domi- 
nate. Originally the neutrino flux formulas given above were derived assuming 
a flat atmosphere, and are therefore only valid up to nadir angles of < 70°. In 
order to describe the full range of zenith angles, 0° — 90°, the cos(6 ) ) depen- 
dence of the it and K critical energies valid at zenith angles below 60° can 
be replaced with a cos(#*) dependence as in [18J. In this work a convenient 
parameterization of this substitution is given: 



cos 2 9 + pi + pi ■ cos P2 9 + p-i ■ cos P4 9 ^ 
1 + Po + Pi + P3 

as given in [TT] and parameters are from the same paper (see eq. [2]). For muon 
neutrinos the low energy function described in eq. [T] and the high energy 
functions have been connected at 500 GeV. The parameters that have been 
tuned based on the Monte Carlo calculation in the region 500 GeV - 100 TeV 
are the normalizations of the pion and kaon terms, respectively A u and B u 
and the spectral index 7. The parameters of the function in eq. [2] are given in 
Tab. [3] and differ for neutrinos and antineutrinos. 



3 Calculation of the muon number of events accompanying neutri- 
nos underground 

The IceCube detector was recently completed and complemented by the Deep- 
Core array dedicated to enhancing IceCube capabilities for neutrino detection 
in the 10-100 GeV energy range. The filter to select internal events vetoed by 
external layer of PMTs is described in [16] ; where the effective area is given. 



COS0* 
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In Ref. [15] the calculation of the efficiency of the veto has been applied focused 
to neutrino astrophysics studies: neutrinos directly produced in sources could 
be distinguished by neutrinos produced in the atmosphere by vetoing events 
where also a muon is detected together with the muon-induced neutrino inside 
DeepCore. We reverse the problem here: we want to keep both the external 
muon and the internal muon-induced neutrino to measure the difference in 
time described in the Introduction. Hence, we calculate the muon event rate 
that can be detected in IceCube and that would accompany a neutrino that 
can induce a muon detectable in DeepCore or in the rest of IceCube provided 
that the track is not starting outside IceCube. We consider pion and kaon 
decays into muon neutrinos and muons and ignore other minor contributions, 
for instance of three-body kaon decays. The kinematics of this two body decay 
is described in Ref.s [T4"fT5] . We assume that the standard kinematics can 
be used for the superluminal neutrinos. This is in general not true, but this 
assumption is justified experimentally from the agreement of the experimental 
data of the two body pion and kaon decays with the standard kinematics. 

In this paper the separation between the muon and the muon induced neutrino 
is calculated and it is about lm (0.1 m) for neutrino energies above 1 TeV 
(10 TeV) for pion primaries and 10 m (1 m) for kaons. These values may 
dependent on the zenith angle. Tracks that are separated by few meters in 
IceCube are a challenge for reconstructions but if the neutrino event may be 
recognized in the inner denser core of IceCue and the muon entrance in the 
detector is detected in the external part the measurement is possible. The 
fraction of primary energy taken by muons and muon neutrinos depends on 
the mass of the primary. The distribution of neutrino energy in the decay 
in the atmosphere are flat over the kinematically allowed regions between 
< E v < Ei(l — 7$), where Ei is the primary pion or kaon energy and 
r i = m ^/ m i = 0.573 for 7T and 0.046 for K ± . Hence the remaining energy 
goes to the muon and it is limited by the condition E^ > E v x [l5j. For 
kaon decays this fraction is much smaller than for pion decays. The fraction 
of neutrinos from kaons and pions as a function of energy is mildly dependent 
on the angle up to 9 = 60°. We took these values from Ref. [19] . 

The minimum muon energy from the decay corresponding to a neutrino energy 
E u has to be enough to be detected underground. The minimum muon energy 
to penetrate at a depth X (in km w. e.) was taken from Ref. [17] and differs 
a little from what quoted in Ref. [15] : E^ min = 0.57 x [(exppf/2.1) — 1] TeV, 
where we take X = h/cosd*, h = 1.95 km w.e. is the depth of the middle point 
in IceCube and cos 9* is given in eq. [3] and accounts for the Earth curvature. 
This expression for the minimum muon energy comes from an average estimate 
of energy losses of muons in ice including stochastic effects from montecarlo: 
- < dE/dX >=a + bEp, where a = 0.268 GeV/m w.e. and b = 0.470 x 10~ 3 
m w.e. -1 . Integrating dX/dE, the inverse f the energy loss, over energy one 
obtains the well known formula E^ min = a/b x (e bX — 1) [14] . 
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We show in Fig. [T]the number of atmospheric muon neutrino events detected in 
IceCube as a function of cos 9 where we consider an energy dependent effective 
area including the effect of DeepCore [TBJ and extrapolate to 86 strings the 
effective area given in Ref. [22] that includes the efficiency for muon detection. 
We neglect dependencies on the angle of the effective area. We show also the 
number of accompanying muons from kaons and pions. The total number of 
events is about 3231 per year. The numbers as a function of cos 8 are given 
in Tab. H] as well. The measurement of the time difference between muons as 
a function of the angle would test the energy dependence of the effect seen 
in OPERA. As a matter of fact the horizontal muons need to have higher 
energies than vertical ones. Since the veto efficiency is about 40% at 1 TeV 
we can expect that in order to recognize the neutrinos these numbers will 
be further reduced but in a way that requires a full simulation for a correct 
estimate. 

For detectors that have good tracking capabilities of multiple muons in a bun- 
dle or for muons in a bundle that have high enough separation, the efficiency 
for this search can be improved not applying any veto requirement to iden- 
tify muon induced neutrinos. The times of arrival of muons in a bundle can 
be measured. Due to the high energies involved, we expect sub-nanosecond 
delays between muons, even if produced by different decays of mesons in the 
shower. The average distance between muons was measured by the MACRO 
experiment [25] at the Gran Sasso laboratory depth and changes from 7 meter 
for vertical incidence up to 19 meters for cosd = 0.55. This average separation 
is small compared to the IceCube dimensions, nonetheless muons with higher 
separations exist and can be detected still with significant statistics given the 
large dimensions (see eg. [23]). 

Very roughly the probability that a muon produced by a neutrino is accompa- 
nied by another muon from another meson that decays in the same shower is 
of the order of 5 — 10% (the typical ratio between double and single muons in 
a muon bundle). This is true if the neutrino energy is of the order of E^ m i n . In 
Ref. |24j it is estimated that the double muons are about 10% (5%) of single 
muons at 3 km w.e. of depth for primary directions of 9 = 0° (6 = 70°). Small 
detectors like MINOS, OPERA (and in the past MACRO) have limited pos- 
sibilities to identify downgoing muons from neutrino interactions. But, if the 
time resolution is very good, the study of the time difference of the events with 
two or more muons could be interesting even without the neutrino signature. 



4 Conclusions 

We have made an approximate calculation of the number of muon events ac- 
companying neutrino events in a detector like IceCube. The determination of 
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Fig. 1. The upper curve is the number of downgoing muon neutrino events per 
year in 86 strings of IceCube as a function of angle, the other two curves are the 
number of accompanying muons from pion and from kaon decays. On the left we 
have assumed the depth of the core of IceCube of 1.95 km w.e. and on the right 
the top of the array depth of 1.5 km we. In the first case the total number of muon 
events is about 3231 and in the second it is 5851 muon events/yr. 
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231,600 


3230.8 


117.9 


3348.7 



Table 4 

Number of atmospheric muon neutrinos [13] and corresponding detectable muons 
per year accompanying a neutrino in IceCube per year at a depth of 1.95 kmwe 
(the core of IceCube). The total in the last row is integrated over the downgoing 
hemisphere. 

the arrival time of the muons and the difference of the arrival of the accom- 
panying neutrinos could be used to verify in the atmospheric neutrino sector 
the OPERA observation. Also the time differences between muons in a bundle 
may indicate that a part of them may have arrived a little earlier. This would 
be the case for the muons induced by neutrinos. 
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